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ABSTRACT 

Experience shows that a vertical lift gate with a flat 
bottom exhibits strong flow induced vibrations and that there 
is an interaction between the vibration of the solid 
boundaries and the exciting force but the nature of these 
vibrations and of the exciting force have not been studied in 
detail, moreover, the interrelationship between exciting 
force, flow pattern, and vibration amplitude have also not 
been set up in general form. 

This thesis, based on extensive experimental work, 
presents the result of an investigation of interrelationship 
of dynamic force and oscillations of a vertical lift gate 
submerged under free surface flow. The primary aim of this 
work is to develop a methodology for separate measurement of 
exciting forces, acceleration, and displacement and to 
correlate the gap ratio with the amplitude of vibrations. The 
main parameters governing flow induced vibrations are gate 
geometry, gap ratio, upstream water .epth, downstream water 
depth, mass, damping, and stiffness characteristics of the 
gate . 


The experiments were conducted in an open water flume 
where the aforesaid parameters were measured with the help of 
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transducers. The signals so measured were amplified by a 
charge amplifier and stored in oscilloscope. The spectral 
analysis of the excitation force and gate response signal was 
conducted using the FFT algorithm. The normalized power 
spectra for the force on a rigid gate and the response of an 
elastic gate throw light on the nature of vibrations of a 
vertical lift gate. The critical vibration were observed near 
gap ratio 0.65 with downstream water depth equal to gate 
th i ckness . 
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LIST OF NOTATIONS 

The following notations are used in this thesis, 
b = viscous damping coefficient 
= critical damping coefficient 
C = linear spring coefficient 
d = gate thickness 
E = modulus of elasticity 
f = dominant vibration frequency 
f = natural frequency (in air) 

F = net fluid force on vibrating gate 
F.^= out-of-phase component of F 

D 

F^‘^= in phase component of F 

F^ = gate mass force component 

F^ = fluid-force on stationary gate 

F^ = resultant fluid force component 

F = Froude number 
r 

g = acceleration due to gravity 

Cj_ = fluid force coefficient 

Cp.^= in-phase force coefficient 

Cp.j^= out-of-phase force coefficient 

C = in-flow added mass coefficient 
m 

hj = upstream water depth 
h^ = downstream water depth 
I = moment of inertia 
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= submergence depth over gate bottom at Vena Contracta 

1 = unsupported length of elastic support 

m = structural mass 

m = structural mass ratio 
r 

Q = discharge in flume 

Re = Reynolds number 

s = gap ratio 

S = power-density function 

S* = normalized power density function 

t = time 

t = dimensionless time 
r 

Tu = turbulence intensity 

V = efflux velocity at Vena Contracta 

V = reduced velocity 
r 

y = gate displacement amplitude 

V = relative gate displacement amplitude 
f3 = damping ratio 

= phase angle between Fp and y 
4> - phase angle between F and y 

w = angular frequency 
a - root mean square value 

V = kinematic viscosity of the fluid. 
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CHAPTER I 


INTRODUCTIOM 


The flow of fluids around structures can cause 
destructive vibrations as well as useful ones. The flow 
induced vibrations have become increasingly important in 
recent years because of its great importance in structural 
design. The objective of designer is to control the 
vibrations when it is objectionable and to enhance the 
vibrations when it is useful. Objectionable vibration in a 
structure may cause eventual failure. 

One of the most dangerous vibrations of hydraulic 
gates with overflow and underflow is flow induced vibrations. 
The gate is a nonlinear system which transforms aperiodic 
input into a different aperiodic output. 

There is a lot of data available on flow induced 
vibrations of hydraulic control gates. However, we can not 
explain the mechanism behind the excitation of self sustained 
gate vibrations because of a great diversity of practical 
gate geometries, boundary conditions and complex nature of 
the dynamic interactions between the flow and the structure. 
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The present thesis is limited to the dynamic 
behavior of planar under flow type gates with a flat 
underside in free surface flow having one degree of freedom 
in the vertical direction. This simple configuration has been 
deliberately chosen in order to eliminate various secondary 
effects due to other boundary and flow conditions that may 
obscure the basic excitation mechanism of a rectangular gate. 

In parallel with Naudascher, E. et. al . . works on 
vortex excited vibrations of an underflow gate in 1986 , the 
present thesis is an experimental study of the excitation and 
elastic response of such a gate that is free to vibrate in 
vertical direction. 

The details are presented in the following chapters 
. Chapter II contains a review of the literature .References 
mainly contain the works of Naudascher, Blevins, Campbell, 
Abelev and Hardwick . In chapter III, a dimensional analysis 
is carried out for the gate vibration problem. In chapter IV 
the description of the experimental setup along with the 
instrumentation used is presented. Results and conclusions 
are presented in chapter V & VI, respectively. Measured 
signals for various gate openings are presented in Appendix - 
A . Appendix - B and Appendix - C contain information on data 
acquisition, signal processing and spectral analysis. 



CHAPTER II 




LITERATURE REVIEW 


2,1 BACKGROUND AND AIMS OF THE STUDY 

Most of the flow i nduced- v i brat i ons can be traced 
to the instability of the flow itself. Whatever be the nature 
of the in stability, its effect in combination with 
disturbances will lead to development of fluctuations of 
velocity and pressure in an initially steady flow, unless the 
Reynolds number of flow is so small that viscous effects are 
large enough to damp this process. Depending upon the 
strength of mechanism by which <i> the intensity of these 
fluctuation get amplified, <ii> their correlation in space 
increases, and (iii) the distribution of their energy becomes 
concentrated around a dominant frequency, an effective force 
fluctuation along the flow boundaries may result and quite 
possibly lead to structural vibrations. 

Muller was the first person who investigated the 
weir gate vibration. Hydraulic gates may be subjected to 
various modes of self excited vibrations differing in their 
physical nature and onset conditions. The mechanism of 
excitation of these vibrations is governed mainly by the 
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character of the nonstationary hydrodynamic forces dependent 
on the motion of the gate on flow. According to Abelev S.A., 
Cl] we assume that this excited vibration of the gates may be 
divided into two basic categories : 

<a) Vibration due to vorticity associated with the 
separation of the flow from the lower lip of the gate when 
the eddy formation in the water past the gate is synchronized 
with and controlled by the structural vibrations ( the eddy 
mechanism of the excitation). 

<b) Vibrations which may occur at high velocity of 
the jet flow oriented along the vertical surface of the gate 
< the jet flow mechanism). 

According to later references CIO] flow induced 
vibrations may be classified into <i> Extraneous ly- Induced 
Excitation, caused by a pulsation in flow or pressure which 
is not an intrinsic part of the vibrating system, <ii) 
Instab i 1 i ty- i nduced excitation, brought about by flow 
instability, and (iii) Movement- Induced Excitation, due to 
exciting forces which arise from the movement of a body 
within the vibrating system. 

When a free shear layer lies close to a stationary 
solid boundary, such that the velocity field of eddies induce 
intermittent reattachment of the separated flow to the 
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boundary. Locher [91 has shown for this case that the 
fluctuating flow pattern exerts a pulsating force with a well 
defined frequency on boundary. 

Kolkman C8], Campbell C4], and Pertrikat C12] 
extended these studies and observed the strong vibrations of 
submerged gates with small openings. They came to the 
conclusion that the vibrations were self controlled and 
originated from the interaction between the unstable shear 
layer separating from the leading edge of the gate bottom and 
the elastic gate. 

Hardwick D.J., [6] studied the vibrations of the 
solid boundaries and their exciting force in detail. In his 
experimental work of excitation and gate response of one 
which is free to vibrate in the vertical direction, he 
measured the dynamic net force acting on the gate bottom and 
the vertical gate displacement. Hardwick established an 
interrelationship between the exciting force, flow pattern 
and vibration amplitude in a nonlinear amplitude response 
curve. However, no spectral information of vibration 
characteristics was given. 

Shyam Sunder C14] studied the gate vibrations using 
a strain gauge bridge for measuring the deflection of the 
elastic cantilever to which the gate was connected. As the 
strain gauge bridge was calibrated by static loading, the 
strain gauge bridge output was proportional to the force 
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acting on the gate as well as displacement of the 
gate . Naudascher , et . al . , Cll] made a comprehensive study of 
the gate vibration problem. Spectral analysis of excitation 
force and the gate response was performed with Fourier 
Analyzer ( HP 5451 C ) for measuring duration of over 20 
seconds. The results were plotted in the form of normalized 
force and the response auto-power spectra. 

The main aim of the present study is to develop a 
methodology for signal medsurement and spectral analysis by 
FFT program on HP main frame computer instead of using a 
Fourier Analyzer. The present experimental study also 
includes the separate measurement of the exciting force, 
vertical displacement of the gate, and acceleration. The 
normalized power density function (with the help of FFT 
program) of the above forces has been found out for both 
elastic and rigid gate case. 


2.2 ORIGIN OF EXCITING FORCE 


In the present study the flow pattern was 
correlated with oscillatory motion of the gate with 
precision. The vorticity of the free shear layer becomes 
concentrated in the growing disturbance as a new eddy beneath 
the upstream third of the gate bottom. 



7 



(a) Maximum upward force (b) Maximum downward 
(gate rising) force (gate falling) 


FIG. 2.1 FLOW PATTERN RELATED TO EXCITING FORCE 


^Higher 
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In the absence of the instantaneous pressure 
distribution on the bottom of the gate some idea of the 

variation of force can be had from consideration of a 
theoretical model of a vortex and its image in a steady flow 
[6]. The unsteady flow has been represented elsewhere by a 

series of steady states, with the growth of circulation of 
each member of the vortex pair equated to the estimated 

circulation in the free shear layer. The dimension of the 
vortex must increase if the vortex is to remain stable in 
flow, and this growth corresponds to the expansion of the 
eddy. The expansion is more rapid as the velocity of the jet 
V^beneath the gate increases. The pressure distribution is 
shown in fig. 2.1. In the fig. 2.1(a) the higher velocity of 
the new eddy gives rise to locally lower pressure on the 
bottom of the gate. Downstream of the new eddy where the 
velocity is lower, the pressure rises to higher value of the 
tail water zone. The maximum upward force is thus applied 
when the new eddy is still small and the old eddy has been 

shed downstream. The gate rises in response to the applied 
force. In fig. 2.1(b) the maximum of downstream force is 
applied when the eddy is large and the zone of reattachment 
is small and the gate starts to fall. By repetition of rising 
and falling motion of the gate, it starts vibrating. 



y/^/'/’/'yyyyyy/’ y / yyyyyyyyyyyyyyyy 


FIG. 5.1 DEFINITION SKETCH 
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CHAPTER - III 

DI MENSI ONAL-ANALYSIS 


The following discussion is based upon similitude 
analysis of Thang and Naudascher [111 The equation for model 
in fig 5.1 is a second order linear ordinary differential 
equation with constant coefficients. 

my+ by + cy = F(t) .... ...<1) 

Here m, b and c denote the structural mass, damping and 
stiffness (in air), respectively, y and its two derivatives y 
and y denote the displacement, velocity, and acceleration of 
vertical gate and F(t) is the fluctuating component of the 
fluid force in the vertical direction acting on the gate 
bottom. 

The above equation may be rewritten in a more convenient form 
by defining, 

c / m = con and b/m = 2Ccon or C = b/2'Ic*m 



u 


where con ( = 2 n fn) is the natural circular frequency of the 
system and C is called the damping factor. Since m, b and c 
are positive C is also a positive number. 

2 

+ (On y = F ( t ) 

+ Y = O V®/ 8 (2) 

r r 

The nond 1 mens 1 onal variables that have been found 
to be useful in describing the vibration of an elastic 
structure in a subcritical (IF < 0.3) steady flow are, 

<i) Geometry (a/b ,s/d ,h /d ,h /d) 

1 2 

(ii) Reduced velocity (V = V / f d> 

r n 

<iii) Damping ratio (/?=b/2(o m) 

r» 

2 

(iv> Mass ratio <m = m/ iodl> 

r 

(v) Fluid force coeff. <C =F/.5pV d^ 1 ) 

F 

(vi> Reynolds number (R« = Vd/v & Vs/r> ) 

The fluid flow and the structure are interactive 
system, and their interaction is dynamic. These systems are 
coupled by the force exerted on the structure by the fluid. 
Because of this coupling no linear solution is possible but 
the gate - displacement amplitude can be expressed in a 
general way as a function of three main groups of parameters 


or , 


y + 2 C y 

Y + 2 /5 Y 


as fol lows s 
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Yo = Function < Boundary geometry. Fluid flow dynamics. 

System dynamics ) 

= Function < s/d , h^/d . , Re . Tu . Vr ,/? , irr > 


It IS difficult to find an analytical expression 
for the relationship between Yo and these parameters. Vickery 
and Watkins assume a linear mechanical system and harmonic 
steady-state force and response function : 

Cf = Cro Sin <ot 

Y = Yo Sin<cot-0) 

CFb V? 

Y = -i- 

Ks 

CFb = Cfo sin 4 > 

Here 4 > is the phase-lag angle. CFb, and Ks denotes out of the 
phase fluid force coefficient and mass damping coefficient 
and {ji>n/co is the ratio of natural frequency (in air) to 

j 

vibration frequency (in water). 


Wn 


a> 


< 5 > 


C*>n 

a> 


[ 


1 + 


Cm 

mr 


] 


Where Cm = m/pd*l is the added mass coefficient. So equation 
(4) may be written as 
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Yo 


1 


An 


Crh Vr* 


Kr 



Cw 


mr 


In addition. Den Hartog, J.P. [53 provides flowing 
formula for natural frequency in still water and coefficient 
of critical damping, 


Wn 



OJn 

fr. = 

Zn 


fn = natural frequency of the gate 
m = mass of the system 
k = spring constant of the system 


5EI 



Q 2 

Here E = 70.0 x 10 N/m for aluminium. 


And critical damping coefficient (be) 


2ncon 
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b 


be 


Where 


1 

c 

Zn 


Xn - Xn+l 


Xn 


> 


Xn 

Xn+l 


n th amplitude during the vibration 
<n+l> th amplitude during the vibration 






GLASS WALLED 
GATE / flume 



FIG- 4-2 SCHEMATIC VIEW OF THE EXPERI MENTAL EQUIPMENT 
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FIG. 4.4 PHOTOGRAPH SHOWING THE EXPERIMENTAL GLASS WALLED FLUME 



CHAPTER - IV 


i ? 


EXPERIMENTAL DETAILS 


4.1 EQUIPMENT 

The experiments were conducted in a water flume 
with 5,0 m. long, 46.0 cm wide and 90,0 cm deep glass 
section. A vertical lift gate fabricated out of aluminium 
sheet 1/16" thick was supported vertically on an aluminium 
cantilever. This type of support simulates the elastic 
support of prototype of the gate. Clearance of 5 mm between 
the gate model and the side walls were found necessary to 
minimize corner eddy effects and to maintain a quasi two 
dimensional flow. A low damping guidance was provided for 
motion in the vertical direction by two roller bearings 
connected to gate model in a vertical plane so that the gate 
had a single vertical degree of freedom. 

The downstream depth in the flume was varied by 
controlling the downstream gate fixed in the flume .■ Di scharge 
was computed by measuring the head of the water over a 
rectangular sharp crested weir in the discharge channel. 





FlG.4-6 BLOCK DIAGRAM OF EXPERIMENTAL PROCESS 







4.2 INSTRUMENTATION 


The experimental setup to determine the vibration 
characteristics of the model is shown in photographs 4.3 to 
4.3 and figures 4.1 & 4. 2. The block diagram of experimental 
setup is also presented in fig 4.6. 

An aluminium cantilever beam of cross section 23 mm 
X 10 mm supported a hollow rectangular box type gate model of 
width 43.0 cm. across and 10.0 cm. along the water flow. An 
inductive d i sp lacement transducer (WIO B&K type) was connected 
to the vertical arm of the cantilever. One accelerometer <B & 
K type 4381 ) was mounted with screw connection on the plate 
provided on the gate. A thin layer of grease was applied on 
the mounting surface before tightening down the accelerometer 
to improve mounting stiffness. The tapped hole in the 
machined part was such that the stud was not forced in to the 
base of the accelerometer. Two piezoelectric force 
transducers (B&K type 8200) were mounted sandwiched between 
two plates on the top of the gate model. The cable used with 
the accelerometer was shielded to prevent noise and clamped 
on the surface The signal from the transducers were amplified 
and stored in a storage oscilloscope (HP type 34301A) . 

With this arrangement it was possible to measure 
the vertical gate displacement Y(t) and the dynamic component 
of resultant vertical force FR(t) acting at location of the 
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F = Fluid-fare 

' o 

F = Resultant -fore 
vibrating gate 


FIG. 4.7 FORCE VECTOR DI 



vibrating gate where the force transducers were mounted, 
When the gate was free to vibrate in the vertical direction. 
From the free body diagram of the gate from below the 
location of force transducers, it would be clear that vector 
sum of the net exciting fluid force F and the internal 
reaction force of vibrating gate mass le , 

Fr = Fo + Fm (with Fm = my ■ ) 

From the knowledge of the Fourier transformed value 
of the amplitude (Yo> . the frequency (w) , and phase angle (^) 
for the dominant vibration component, it was possible to 
determine the magnitude Fo and the phase angle 4 > of the net 
fluid force component of the same frequency to by method of 
the force vector diagram (fig 4.7). 

4.3 PROCEDURE 

A gate model of d = 10.0 cm. and weight 8.25 Kg was 
used . The gate was connected to a support by means of an L - 
shaped cantilever. A number of holes were made on the 
vertical arm of the flat to enable setting of the desired 
gate opening. A cantilever span of 15 cm was used. For each 
gate opening, the downstream depth was varied with the help 
of a downstream control gate. The discharge was measured with 
the help of a sharp crested full width weir. The velocity of 
the jet Vj underneath the gate was calculated as follows : 
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1 . 45 

Q = 0.857 h m®/ sec 

where h = head over the crest of the rectangular 

weir (m) and 

Vj = Q / area of gate opening 

Fig 4.2 shows the schematic view of the 
experimental equipment. Various photographs can also be seen 
to have a clear picture of the setup. 

For signal analysis the transducers were connected 
to charge amplifiers for amplification of signals and the 
amplifiers were connected to a storage oscilloscope. This in 
turn was connected to a PC in which a GPIB interface card 
was installed. The signal was first stored in the storage 
oscilloscope and then transferred to the PC. 

For determination of damping characteristics the 
displacement signal was stored for various depths of still 
water in the flume. The deflection of the model gate under 
various loads was found out with the help of the displacement 
transducer . 


Spectral analysis of the resultant force and gate 
displacement signal was performed with FFT software on the HP 
main frame computer. The results were plotted in the form of 
normalized force and displacement auto-power spectra. The 
phase angle <p between the resultant force and the gate 
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displacement was obtained 
software . 

4» 4 SPECI FI CATI ONS 

4.41 OSCILLOSCOPE 

Model number : 

Manufacture : 

No. of input channels : 

I nput coupl i ng : 

Storage mode : 

Normal mode : 

Vertical sensitivity : 

Vertical gain accuracy (dc> ; 
Time base accuracy j 

Time base range : 

Maximum sample rate : 

Memory depth : 

Signal averaging : 


from cross-spectral analysis 


HP 54501 A 
Digitizing oscilloscope 
Hewlwett Packard 
4 

AC . DC 

Oscilloscope digitizes the 
input signal before 

displaying the trace 

Behaves like a conventional 
analog instrument 
Maximum 5 mv/div 
Minimum 5 v/div 
+ 1 . 5 « 

0.005 % 

2 ns/div to 5 s/div 
10 Msa/s 

501 P.oints (display) 

1024 Points (via HP-IB) 

1 to 1024 
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4.42 AMPLIFIER FOR DISPLACEMENT TRANSDUCER 


Type : 
Manufacturer : 
Measuring range : 
Measuring rate : 
Input & Output control signals : 
Linearity division : 
Interfaces : 

4.43 CONDITIONING AMPLIFIER 


ALPHA 3000 

Hottinger Baldwin Messtechnic GMBH 
± 10 V 
25 Hz 

T T L level 

< 0.01 % 

IEEE 486-78 RS 232-C 


Type 

Manufacturer 
Charge input 

Maximum input 
Sensitivity conditioning 
Amplifier sensitivity 
Signal Output 

Maximum Output 
Frequency range 


Charge Amplifier Type 2635 
Bruel & Kjaer ; Denmark 
Via 10 - 32 NF and BNC coaxial 
socket 
ac 10®. pC 

0.1 to 10.99 pC / ms"* 

0.01 mv to 10 V / pC 

Via 10 - 32 NF and BNC coaxial 
socket 

8 V < 8 mA > 

2 Hss to 100 KHz (Acceleration) 

10 Hz to 10 KHz (velocity) 

10 Hz to 1 KHz (Displacement) 
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Power requirement 


4.34 ACCELEROMETER 


: Three 1.5 Alkaline cells or +6 to 
+28 V single or ± 14 dual polarity DC 


Type 

Manufacturer 
Charge sensitivity 
Voltage sensitivity 
Frequency range 
We i ght 

4.35 FORCE TRANSDUCER 


B & K type 4581 
Bruel & Kjaer ; Denmark 
: 10.0±2.0«pC/ ms"^ 

: ia 8 . 0 mV / ms * 

4.8 KHz 
43.0 grams 


Type 

Manufacture 
Force range 
L i near i ty 

Charge sensitivity 
St i f f ness 
We i ght 


: B & K Type 8200 

: Bruel & Kjaer ; Denmark 

: 1 KN tensile to 5 KN compressive 

: < ± 1.0 % of maximum force 

: 4.0 pC / N 

; 5.0 X 10® N / m 

: 21 gram 


4.36 INDUCTIVE DISPLACEMENT TRANSDUCER 


Type : W 10 K 

Manufacturer ; Hettinger Baldwin Messtechnic 

GMBH 


Nominal displacement 


± 10 mm 
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Nominal sensitivity 
Nominal Output signal span 
Linearity Deviation - 
Carrier frequency 
Maximum cable length 
Weight (with plunger) 


± 80 mv / V 
160 mv / V 
< ± 0.4 % 

5 KHz 
900 m 
40.0 gms . 
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CHAPTER -V 


RESULTS AND DISCUSSION 


The results are tabulated in Table 1 to 3. The 
recorded measured signals are shown in Appendix - A . The 
spring constant of elastic support is calculated both, by 
equation given in the chapter III and by load deflection 
curve. It has come out as 117.64 x 10^ N/m. The coefficient 
of critical damping and coefficient of damping for various 
depths of still water from bottom of the channel is given in 
Table ^1. It is observed from Table // 1 that the coefficient 
of damping decreases with increase in depth of still water in 
the flume. 


The influence of various parameters on vibration is 
shown in Table if Z. The Reynolds number s / v lies between 
143000 to 180000 and one interesting thing observed here is 
that this Reynolds number increases with increasing gap ratio 
(s/d) but decreases with increasing downstream water depth 
for the same gap ratio. 

Table if 3 throws some light on the relationship 
between the gap ratio and the mean amplitude of the measured 
signals of displacement and force. This mean value of gate 



displacement and excited force is averaged over a sufficient 
time interval. The influence of gate opening on force and 
displacement magnitude can be seen in Table if 3. 

Results for the phase angle between fluid force and 
gate displacement are presented in fig. 5.5. As the gap ratio 
increases the phase angle <!> is observed to decrease starting 
from 90 degree. In order to find the most critical flow 
situation with respect to vertical gate vibrations the 
normalized power density function of displacement vs 
frequency is plotted in fig. 5.6. It is observed with the 
help of fig. 5.6 that the peak of the spectra for the free 
vibrations match with the peak of spectra for gate vibration 
with gap ratio 0.65 . Also from Table # 3 we see that the 
maximum mean displacement occurs for the case of gap ratio 
0.65 Hence we can say that the most dangerous vibrations were 
found at the gap ratio 0.65 . 

For understanding the nature of vibrations, the 
power spectra of the force on a elastic gate <Sff> and the 
response of an elastic gate (SYY)are plotted in fig. 5.7 to 
fig5.12. These spectra are plotted in a normalized form using 
normal scales for various gap ratios and downstream water 
depth. The linear dependence of dominant excitation frequency 
on flow velocity and gap ratio for the rigid gate and the 
resonance character of the response for the elastic gate are 
clear indication of vortex excitation . 
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For the clear visualization of the vortex pattern 
underneath the gate for various gap ratios and downstream 
water depths photographs have been taken and shown in fig. 
5-13. Visualization of the shear flow zone underneath the 
gate reveals that the broad-banded peaks in these force 
(elastic gate) spectra are associated with large scale 
vort ices. 
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CHAPTER -VI 

CONCLUSIONS 

The vibration characteristics of a vertical lift 
gate used to control flow depend on the geometry of the gate, 
upstream and downstream water depths, gap ratio, mass of the 
gate and damping characteristics of the gate. Using 
dimensional analysis, a functional relationship between the 
dominant amplitude of vibration and the above listed 
variables has been obtained. The new methodology for separate 
measurement of force, acceleration and displacement with help 
of transducers, amplifier and oscilloscope is used 

successfully. Spectral analysis of measured signals with the 
help of FFT software provides good knowledge about the nature 
of vibrations of vertical lift gate. The following 

conclusions can be drawn from the experimental study : 

(i) Maximum vibration of the vertical lift gate occurs 
at a gap ratio of 0.65. 

(ii> Vibration amplitude reaches a maximum at a 
submergence over the gate bottom ratio of about h^ /d - 1.0 
and decreases continuously with increase in the h^ /d ratio. 
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(ill) The nature of the vibration can be explained with 
help of the power spectra of force on a rigid gate (S > and 

FF 

of the response of an elastic gate (S ). 

VY 

(iv) As the velocity increases beyond the vortex 
resonance range more and more energy is transferred from the 
fluid to the body. 

The above conclusions are similar to those of 
Naudascher , E . , et . al . , [113 . 

6.1 DESIGN CRITERIA DERIVED FROM THE EXPERIMENT 


Based on the conclusions the following design 
criteria can be formulated : 

<i) Resonance freq,uencies of gate hoist systems and gate 
members should be well above the dominant frequencies. 

(ii> Unstable flow separation and / or reattachment zone 
should be prevented. 

(iii) Water surfaces unstable because of turbulence, 
waves, hydraulic jumps, water column separation, should have 
no opportunity to hit the structural surface. 
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APPENDIX - B 


DATA - ACQUISITION 


Data Acquisition generally relates to the process 
of collecting the desired process in digital form, as 
rapidly, accurately, completely and economically as possible. 
Some form of storage is generally required to permit 
subsequent processing and, indeed, to carry out reprocessing 
at later time when different analysis procedures are seen to 
be required. For this purpose a data acquisition system which 
acts both as a buffer and a converter between the data source 
and analysis operation is required. 

The real world events are continuous (or analog) in 
nature but computer is a digital device. Hence it requires a 
conversion of analog signal to digital form samples are 
quantified for digital processing by converting the sampled 
value into a equivalent numerical value by an analog to 
digital converter, each value is represented by a finite 
number of ON / OFF states of the solid state circuit. 

In a data acquisition system the following 
operations are involved - generation of input signals by 
transducer, signal conditioning, data conversion < from 
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analog to digital and vice versa ) data storage and display 


and finally. 

data 

process i ng . 

In 

the 

s i gnal 

condi t ioner 

operation we 

have 

to amplify 

the 

weak 

i nput 

signal from 


transducer with the help of a charge amplifier. In this 
process we also filter out the unwanted frequencies 
originated from the transducer .The analog transducer output 
is converted into digital form using a analog to digital 
converter (A/D) built into the measurement system before data 
transmission to a personal computer. 

After A/D conversion the digital data is stored in 
the hard disk of the PC and subsequently in floppies. 
Standard interface are amicable for controlling digital 
instrument and transferring data to the PC. We have used 
general purpose interface bus (GPIB) for this purpose . The 
bus protocol was evolved by Hewlett - Packard company (USA) 
under the name of HP-IB for the interconnection of 
instruments . The GP IB interface card provides a low cost 
interconnection for up to 15 digital instruments over a 
distance (about 20 mts). The GPIB consists of 16 lines of 
which 8 are meant for parallel data transfer and 8 form the 
command channels. The maximum data transfer rate is about 1 
MByte / sec, but transfer rates achieved in PC based 
instrument is about 250 MByte/sec . A GPIB cable is 
terminated at both ends by identical connectors which permit 
Piggy - back connections . 
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In the present work, a GPIB interface incorporated 
with HP digital storage oscilloscope is used for transferring 
of data from the oscilloscope memory to computer. For IBM-PC 
and a compatible GPIB interface called GPIB-PC2 is used. It 
is possible to use this software through user programs 
written in FORTRAN and PASCAL. Commands for sending data from 
the HP oscilloscope to PC-XT using the IBIC menu are given 
below : 


IBIC 

1 IBFIND Osc 

2 IBWRT ACQUIRE: TYPE NORMAL; COUNT 1; POINTS 512” 

3 IBWRT DIGITIZE CHANNEL 1" 

4 IBWRT SYSTEM: HEADER OFF;:EOI ON” 

5 IBWRT WAVEFORM: SOURCE CHANNELl ; FORMAT ASCII” 

6 IBWRT " :WAVEFORM:DATA?” 

7 IBRDF OUTPUT 

IBCLR 

I BLOC 
QUIT 


The first command is to open the oscilloscope. The 

second command is the acquire subsystem. Third and fourth are 

» 

for counting 512 data and its digitization . Command six for 


format for source and output 

. Seventh 

command 

br 1 ngs 

the 

osc i 1 loscope 

in action to get 

ready to 

transfer 

data . 

The 

next command 

stores the data 

in PC under the 

file 

name 


OUTPUT. 
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C 


C 1 FAST FOURIER TRANSFORMATION 


Fast fourier transform is a very efficient 
algorithm which calculates the discrete fourier transform 
(DFT) based on discrete version of the fourier integral pair. 
From this it is possible to make rapid transformations 
between time and frequency domains. 

The DFT is a finite sequences f Xr } , r = 

0,1, 2, 5 <M-1) is a new finite sequence C Xk > 

def 1 ned as , 

, i - .t / N ) 

" r=o ■■ K = 0,1 .2 <N-1) 

r 


If we are working out values of Xj^ by a direct 
approach we will have to make N multiplications of the form C 
X } X for each of N values of X, and so the 

r k 

total work of calculating the full sequence Xj^ would require 
N* multiplications. The FFT reduces this work to a number of 
operations of the order NLog^N. The FFT therefore offers an 
enormous reduction in computer processing time and also 
increases accuracy. 
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The signal to be analyzed is initially in analog 
form and is (i) sampled <ii> digitized (iii) stored in the 
computer memory and <iv) spectral analysis is carried out 
using fast fourier transforms <FFT> . 

C 2 SPECTRAL ANALYSIS 


Spectra provide a frequency-domain resolution of 
physical quantities such as displacement. velocity, 
acceleration and force. Spectral data can be processed in 
such a same manner as time domain data in determining 
correlation functions. There are two classes for spectral 
data : 

<i> CROSS-SPECTRUM 

The cross spectrum measure the similarity between 
two signals in frequency domain, it is the fourier transform 
of the cross-correlation and expresses the similarity as a 
function of frequency. Conjugate multiplication is applied to 
calculate the cross spectral density function. 

(ii> AUTO-SPECTRUM (POWER -SPECTRUM) 

It is a similar measurement in frequency domain, 
the measured quantity being multiplied by itself such that 



the coefficients of the components are expressed as the 
square of their magnitude. 

An electrical analogy by Betts considers a periodic 
function F<t) which defines either voltage or current 
variation associated with a circuit in which the power is 
dissipated by a one-ohm resistor. This would be regarded as 
normalized power ( based on unit resistance > and thus the 
average value is given by : 

1 +T/2 

P = S f*(t) dt 

T -T/2 

In terms of the fourier transform this can be rewritten as. 


P 


1 +T/2 +T/2 

/ f<t) dt / F(jw> expCjot) df 

T -T/2 -T/2 


P 


1 


T 


+T/2 

S |F<jw>p df 
-T/2 


Since Hz = (1/time) = <l/sec), the integral reduces to 
(watts/Hz) . 


An alternative but complementary way of describing 
the power spectrum is as follows, physically this is simply a 
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measure of how much energy is contained within the signal in 
each frequency band. If X<t.co>,/Nw> is that part of the signal 
X<t> which lies in the frequency band c:. <ji center at 
frequency w. Then the mean square value within this band is 


ioi) 


= 1 i m C- 
T-* 00 


1 ^ 

T 



t ,<o, ^ (o) dO 


Then power spectral density can be defined as 


|F(w)j^ = lim [ S^<w)/( o> X >] 

f-* 00 

where x^is the mean square <or variance o'* ) of the original 
signal X(t) . 
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C PROGRAM FOR SCALING THE SIGNAL DATA 

integer*^ s 

dimension e ( 1 0£4 ) , b ( 1 0E4 ) 
integer^A amax,amin 

open<uniT = 1£. filer's . in' .. sTatus='old' ) 
read(1S,»-){a{ i ), i^l ,51£) 
closetunlt=l2, status*' Keep' J 
emax=0 

amin=l 000000 
do ^ 0 1*1 ,51£ 
aniax*maxO c amax . a ( i > ) 
emln*minO(emln/a( 1 J ) 

1 0 cont inue 

print*,amax,amln 
do £0 i*1 ,51£ 
c=256*< a1 1 )-atnln) 

1 )=c/<amax-amin) 

EO continue 

do 30 1=1,512 
a( 1 )*&l 1 ) 

30 continue 

open t unit *E0 , r 1 1 e* ’ int eg . dat ' , status® 'unKnown • ) 
write{£0, » ) (a( i ) , i®1 ,5ie) 

Close ( uni t=E0 , St at us® ' Keep ' ) 

stop 

end 



PROGRAM ELASTIC 


CALCULATION OF FREQUENCY SPECTRUM OF DIGITIZED 
SIGNALS USING FFT FOR ELASTIC GATE. 

dimension e{3£00), power! 3£00 ) / «ut o ( 3£0 0 )>alag(3£0D) 
dimension b< 3£00 ) , a£( 3£00 ) 
complex x(3£00) 
common / tr /troom 

print *, 'whet are nu ndata del ' 
reed *i nu, ndata, del 


n=£**nu 

del = l 0, 0*del/T 0£4 
f me x» 1 . 0/de 1 
tt=dei*(n-l ) 
pl=3, 14159H6 

sea le®1 . 0 

open(unit«=£l ,file»'integ.det' , status* 'old' ) 
reed<£1,*) ( a( i >, i=1 , ndata ) 
close {unlt=£1 , status® ' keep ' ) 

print *, 'signal data has been read' 

remove stray mean. 

call trap < a , ndat a , e ) 
do 65 i = 1 , ndat a 
aC 1 )«a( 1 )-6 
cont 1 nue 

filter data (Tukey- Henning formula). 

t i 1 *0 . 5»de 1* < ndata- 1 ) 
do 75 1®1 , ndata 
t l = del* < 1-1 j 

wt = 0 . 5* < 1 . 0+cos{ pi* (t 1-tt 1 )/tt 1 ) ) 

aE< i )=a< i )**£ 

e( 1 ) = e! i )*wt* scale 

cont i nue 

call trap ( a£, ndata , var iance ) 

do 45 i=ndats+1,n 
at 1 ) = 0 . 0 
cont 1 nue 



1 00 


1 1 1 


lib 


iE30 

c 

1200 

1250 

c 

1275 


5 


6 


call f f tr ( a , b , X , nu , del ) 

do 100 1=1 ,n 
r=re*l ( x( i ) ) 
s = iinag( x C 1 ) ) 

powerCl )=sqrt(r*»£+s'»:%£) 
cent i nue 

print 'frequency spectrum calculated' 


power 1 =-1.0 E+06 
Do ill i=l , n 

power I =max 0 1 i f 1 X ( power l ) , in x (power! l ) j > 


cont inue 
do 116 i=l,n 
power ci)“powerCi) /power l 
cont inue 


nd2=n/2 
f nyq=0 . 5*fmex 


open ( un it =£4 , f i le= ' f req . dat ' , status* ' unknown ' ) 
write(£4, 1230) 

format C /Ex j ' freq in Hz.; amplitude dimensionless',//) 

writefe^, 1200 ) 

f ormet ( 6x , ' f req ' , 7x , 'ampl 1 tude (real, Imaginary mod)',//) 
write(£'4,l£50) (b(i),power(i),i = l,nd£) 
format (£x ,c1£.5,4x,e1£.5) 
wr ite(E4, 1275) fnyq 

format { /5x , 'Nyquist frequency (Hz) ',e1£.E> 
closeCunit=£4,status='keep' ) 

stop 
end • 

subroutine f f t r ( a , b, x , nu, del ) 
dimension x(3200) ,b(3E00) ,a(3200) 
complex x,u,w,t 
n=2‘t:»nu 

pl=3 . 141592653 
tt=(n-1 )’t=del 

do 5 3 j = 1 , n 
b(jj) = (jj-1)/tt 
cont inue 

do 6 1=1 , n 

x(i )=cmplx(a(l3,0.0) 

continue 

do 20 1 =1 , nu 

{ nu + 1 -1 ) 
le1«le/2 
u*! 1 . , 0 . ) 

w=cmDlx (cos(pl/f loat ( lei ) ) , -sin(pl/f loat ( lei ) > ) 


do EO , le 1 


do 10 iej,n,le 
ip*i+lel 
t = x C i ) + >: ( ip ) 
x(ip)«(x<i)-xCip>)>»-u 
10 x(i)-t 

£0 u=u*w 

nvE=n/£ 
nml =ri-1 

do 30 i =1 , nmt 
i f { i . ge . j ) go to £E 
tax ( j > 

X ( j ) a X ( i ) 

X < i )at 

£5 kanv£ 

35 i f C k . ge . J )go to 30 

jaj-k 
kak/8 
go to 35 


30 j®j+k 

ret urn 
end 

subroutine simpson (x,n,e) 
dimension x(3£00> 
f xao . 0 

do 5 ia£, n-1 ,£ 

5 fx=rx+x(l) 

f yao . 0 

do 10 ia3,n-£,£ 

1 0 f yafy + x { i ) 

e=(x( 1 ) + x(n)-»-4. 0»=fx+e.0»fy)/i3.0*{n-1 ) ) 

return 

end 

subroutine trap<x,n,e) 
dimension x(3£00) 

f x = 0. 0 

do 10 i*£,n-1 
f x=f x+x ( i ) 

10 continue 

ea(xf1) + x(n)+£.0’>'fx)/{£.0»(n-1)) 

return 

end 



PROGRAM RIGID 


C CALCULATION OF FREQUENCY SPECTRUM OF DIGITIZED 

C SIGNALS USING FFT FOR RIGID GATE. 

dimension acc( 1 SOO ) , power ( 1 SOO ) , aut o< ieo0),alag(1£00) 
dimension bi1200),a£(1200), f orce( 1 £00 ) 
complex eccf ( 1 EOO ) , f orcef ( 1 £00 3 

print «= .. 'what is the mass of the gate in kg' 
print ' ' 

read *, am 
print * , ' ' 

print if, 'give the force and acceleration scale factors' 

print *, 'in units of N/mv and m/E*s»mv' 

print*, ' ' 

read *, sfr,sac 

print ' ' 

print * , 'give the oscilloscope voltage scales for the force' 

print »■ , 'and acceleration signals in mv/division' 

print ' ' 

read », scala,scalf 

print * , ' ' 

print *, 'IMPORTANT: THE TIME SCALES ON THE OSCILLOSCOPE FOR' 
print *, 'THE TWO SIGNALS MUST BE EQUAL (IN SECONDS)' 
print *, 'what is the time scale used on the oscilloscope' 
print « , ' ' 

read *, del 

nu=9 
n«=£**nu 
ndat B=n 

del =1 0 . 0*de 1/n 
f max® 3 . 0/del 
tt=del* (n-1 ) 
pi»3. 34159S653 

opentuniteSI ,file®'acc.dat ' ,status='old') 
read(£l,*) (acc{ i ) , i=l , ndata ) 
close (unit =£1 ,statUB='keep' 3 

open ( un it =££, f i I e= ' f orce . dat ' , status® ' old ' 3 
read(££,*) ( force ( i ), i=1 , ndata ) 
clOEe<unlt=£E,statUB='keep' 3 

print »•, 'signal data has been read' 

c filter data (Tukey- Hanning formula). 

t t 1 *0 . 5*del ♦ ( ndat a-1 ) 
do 75 i=1,ndata 
t l=del*( i-1 3 

wt=0 . 5* ( 1 . 0+cos( pi*( t 1-tt 1 )/tt 1 ) ) 

acc( i )=8. 0*(acc( i )/3£640 )*wt*scel a*Bac*am 

force(i 3=8. 0*( force( i )/3£640 ) *wt ♦seal f*sfr 


cont i nue 



65 

c 

45 

1 00 

£26 

1 16 

1 23 0 
1 200 
1250 
1£75 


remove stray mean. 

call trap(acc,ndata,ea) 
call trap<force,ndata,ef ) 
do 65 i=1 iHdata 
acc(i )«acc(l )-ea 
force(i)=force(i)-ef 
cont i nue 

add zeroes to complete string. 

do 45 i=ndata+l , n 
acc ( i )®0 . 0 
force ( i ) *0 . 0 
cont 1 nue 

call f f t r ( acc , b , accf , nu, del ) 
call f ftr ( force, b, forcef, nu, del ) 

do 100 i=1 , n 

r=real(accf(l) -f orcef ( i > ) 
s=imag(accf( 1 )-forcef ( i ) J 
power( i ) = Eqrt <r*=»£+E**£) 
cont i nue 

print **■ , 'frequency spectrum of fluid force calculated' 

powerl *-1.0 E+06 
Do ££6 i *1 , n 

powerl » maxOl if lx (powerl ) , if i X <power( 1 ) ) ) 
cont i nue 
do 116 i * 1 , n 

power { i ) * { power ( l ) )/■( power! ) 
cont i nue 

ndE=n/£ 
f nyq=0 . 5*fmax 

ope nC unit *£4 , f i 1 e* ' f req . dat ' , status*' unknown ' ) 
write{£4, (£30) 

f orniat ( /Ex , ' f req in Hz.; amplitude dimensionless',//) 
write(E4, lEOO) 

f ormat ( 6x , ' freq ' , 7x , ' amp! Itude (real, imaginary mod)',// 
writeC£4, 1£50 ) (bC i ), powerl i ) , i*l ,nd£) 
format (£x,e1£.E,4x,e1£.5) 
wrlte(24, 1275) fnyq 

f ormat < /5x , ' Nyqu ist frequency (Hz) ',elE.5) 
close(unit*£4,status='keep' ) 

stop 

end 

subroutine f f tr ( a , b , x , nu , del ) 
dlmens 1 on x( iEO0),b(l2OO),a(!EOO) 
complex x,u,w,t 


n*£**nu 

nf *-x 1 



do 5 j j * 1 1 n 
b(jj) « (jj-D/tt 
cont i nuo 


do 6 i=1 , n 
X < 1 j = cffipl X ( a ( i ) , 0 . 0 } 
cont inuo 

do £0 1 *1 , nu 
le=E=** ( nu+ 1 -1 ) 
le!=le/E 
u={ 1 . , 0 . ) 

w=cinplx (COS (pi /float ( lei ) ) ,-sln( pi /float (let ) ) ) 
do £0 j=1 , lei 


do 10 i- j , n, le 
ip=i+le 1 
t=x ( i )+x( ip) 

X t ip)® ( x{ i )-x ( ip ) )*u 
X < i )=t 
U = u s*'W 


nv£*n/£ 
nmi ®n-i 
j«i 

do 30 i=1 , nml 

if ( 1 , ge . 3 )go to E5 

t®x ( j ) 

X ( j ) = x ( 1 ) 

X { i )«t 
l5=nvS 

If ( k . ge . 3 )go to 30 

k = k/£ 
go to 35 

3 = 3+1^ 


return 

end 



subroutine simpeon <x,n,e) 
dimension x<1£00) 
f XB'0 . 0 

do 5 l«=E< n-l , E 
5 rx»rx+xci) 

ry=o . 0 

do 10 i=3,n-£,e 
1 0 f yssf y + x ( 1 j 

€=( x( 1 )+x <n)+4. 0*fx+E. 0*fy)/(3. 0*(n-1 ) ) 

return 

end 

subroutine trap{x,n,e) 
dimension x(ISOO) 

f x*0 . 0 

do 10 l=E,n-1 
rx = rx + x c 1 ) 

10 continue 

e=<x( 1 )+x(n)+S,0*fx)/(£. 0*(n-1 )> 

return 

end 
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PROGRAM PHASE 


CALCULATION OF PHASE ANGLE BETWEEN THE FLUID FORCE 
AND THE GATE DISPLACEMENT. 


dimension acc(IEDO), power ( 1£00)/BUto(t£DO}(aIao( tBOO) 
dimension b < 1 £00 ) , a£ (1 £00 ) , force ( 1 £00 ) , disp( ! £00 ) 
complex accf ( 1£00) ,forcef < 1£00) ,dispf < 1£00 > 

print , 'what is the mass of the gate in kg' 
print ' ' 

read am 

print *, ' ' 

print 'give the force and acceleration scale factors' 

print 'in units of N/mv and m/ses*mv' 

print ' ‘ 

read sfr,sac 

print ' ' 

print 'give the oscilloscope voltage scales for the force' 

print =» , 'and acceleration signals in mv/division' 

print ■* / ' ' 

read *, ecala,scalf 

print *, ' ' 

print *, 'IMPORTANT: THE TIME SCALES ON THE OSCILLOSCOPE FOR' 
print 'THE THREE SIGNALS MUST BE EQUAL' 

print *, 'what is the time scale used on the oscilloscope' 
print ' ' 
read del 

nu*9 

n-Eei'nu 

ndata«n 

del=l 0 . 0«-del/n 
fmax»l . 0/del 
tt®dslt=ln-l ) 
pl=3. 1 

open! unit «=£1 , f i 1 e* ' acc . dat ' , status* 'old' ) 
reedC£l,*) ( acc( i ) , i*l /ndata ) 
close{unlt*El , status® ' keep ' ) 

open(unit=£E, file® 'force.dat ' , status® 'old' ) 
reed(££,*) ( f orce( i ) , i®1 , ndata) 
close(unit=E£, status®' keep' ) 

opent unit ®£3, fi le*' disp.dat ' , status® 'old' ) 
re8d<£3,*) ( disp < i ) # i= 1 , ndat a ) 

clOBe<unit®£3,status®' keep' ) 
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print », 'signal data has bean read' 

c filter data (Tukey- Hanning formula). 

ttl «0 .St'del^Cndata-l ) 
do 75 1«1 , ndata 
t l=dei* ( i~l ) 

wt*0 . S* ( 1 . 0+cos( pi*( t 1-tt 1 J/t t 1 J ) 
acc{ i ) = 8. 0*(acc( i )/3S&40>»wt*scsla*sdc*ain 
force ( i ) = S . 0* ( force ( 1 )/3S640 ) *wt *scalf *Brr 
dlsp( 1 }=8. 0*(dlsp( 1 )/32640 )*wt*scalf*sfr 
75 continue 

c remove stray mean. 


call trapiacc, ndata, ea) 
call trapi force, ndata, ef) 
call trap(disp, ndata, ed) 
do 65 i=1 , ndata 
acc( i )*acc( i l-ea 
force ( i ) “force < i )-ef 
dispC i )=di-3p( i )--ed 

65 continue 

c add zeroes to complete string. 

do 45 i«ndata+1,n 
acc( i )*0 . 0 
f orcei i )*0 . 0 
disp(i)“0.6 
45 continue 

call ff lr(acc,b,accf ,nu,del) 
call fftr(force,b,forcef ,nu,del) 
call ff trCdisp,b,dispf , nu,del ) 

do 100 i “1 , n 

r«real ( accf ( i > -force f( i ) 1 
s“imag( accf ( i )-f orcef (ill 
power t i )*sqrt ( r*eS+B**£ ) 

rd^realtdispfli)) 
sd=iiBag{ diBpf ( i > ) 

p1»atan(r/&) 
p£®atan( rd/sd) 


alag( i )*ab-3{p1 -pS) 



1 00 


1 16 


1 E30 
1 £00 
1£S0 
1 £75 


5 

6 


cent i nue 

print », 'frequency epectrum of fluid force calculated' 

print ' ' 

print 'phase angle between fluid force and gate' 

print 'displacement calculated' 

powerlepower( 1 ) 

do 116 i =1 , n 

power ( 1 )* power ( i ) /power l 

continue 

ndE=n/£ 
f nyq=0 . 5* fmax 

open( un it aE4 , file* 'freq.dat ' , status*' unknown ' ) 
write{£4, 1£30 ) 

f ormat ( /Ex , ' freq in Hz.; amplitude dimensionless',//) 
write(£4, t£00) 

formattex, 'freq',7x, ' ampl 1 t ude ' , 8x , 'phase angle ',// ) 
write(e4,1£50) (b(i),power(i),alag(i>,i*1,nd£) 
forinat{£x,eie.5,4x,e1£.5,4x,e1£.5) 
writet£4, 1E75) fnyq 

format < /5x Nyquist frequency (Hz) ',e1£.5) 
close(unit*£4,Etatus*' keep' ) 

stop 

end 

subroutine fftr( a, b, x , nu, del ) 
dimension x(lEOO),b(l£00),a(l£00) 
complex X , u , w , t 


n*£*t nu 

pi*3. 141S9E653 
tt*{n-1 )*d6l 

do 5 j j *1 ( n 

b( j j ) * ( j j-1 )/tt 

cont i nue 

do 6 i*1 , n 
X ( i )=cmpl X ( a( i ) , 0 . 0 ) 
cont i nue 

do £0 1*1, nu 
1 e*E** ( nu+l -1 J 
lei *le/E 
u* ( 1 . , 0 . ) 

w*cmplxt cos (pi/float (lei )),-sin(pl/float(lel ))) 



/b 


do £0 j*! , lei 

do 10 i* j , n, le 
ip*i+Ie1 
t“X{i)+x(ip) 
x{ip)®{x(i)-x<ip))»-u 
10 x(i)=T 

SO u*=u*w 

nv£=n/£ 
nml ^n-1 
j = l 

do 30 i=1 , nml 

if ( 1 . ge . J )go to S5 

t=x { j ) 

X ( j ) = x ( i ) 

X C i ) = t 

S5 k=nvS 

35 in k -ge . J )go to 30 

j = j-k 
k=k/£ 
go to 35 
30 

ret urn 
end 

subroutine simpson (x,n,e) 
dimension x<1£00) 
f x*=0 . O 

do 5 l=E,n-1,E 
5 rx=rx+x(i) 

f y=0. 0 

do 10 i=3,n-E,£ 

10 ry»ry+x{i) 

6=( x{ 1 ) + x(n)+4.0*fx+£.0»fy)/(3.0*(n-1 ) ) 

return 

end 

subroutine trap(x,n,e) 
dimension x(1£00) 

f x = 0. 0 

do 10 l*E,n-1 

f x=f X +x ( i 1 
10 continue 

e«(x(1>+xln)+E.0*fx)/(£.0»(n-1)) 


return 

end 



